Power-Law Decay of Standing Waves on the Surface of Topological Insulators 
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We propose a general theory on the standing waves (quasiparticle interference pattern) caused 
by the scattering of surface states off step edges in topological insulators, in which the extremal 
points on the constant energy contour of surface band play the dominant role. Experimentally we 
image the interference patterns on both Bi2Te3 and Bi2Se3 films by measuring the local density of 
states using a scanning tunneling microscope. The observed decay indices of the standing waves 
agree excellently with the theoretical prediction: In Bi2Se3, only a single decay index of —3/2 exists; 
while in Bi2Te3 with strongly warped surface band, it varies from —3/2 to —1/2 and finally to — 1 
as the energy increases. The —1/2 decay indicates that the suppression of backscattering due to 
time-reversal symmetry does not necessarily lead to a spatial decay rate faster than that in the 
conventional two-dimensional electron system. Our formalism can also explain the characteristic 
scattering wave vectors of the standing wave caused by non- magnetic impurities on Bi2Te3. 

PACS numbers: 73.20.-r 68.37.Ef 73.43.Cd 72.10.Fk 



The discovery of topological insulators (TIs) has at- 
tracted a great deal of attention [THU]. The three- 
dimensional TIs are characterized by the gapped bulk 
states and gapless surface states (SSs), which are pro- 
tected by time-reversal symmetry (TRS) and consist of 
an odd number of spin-helical Dirac cones. Exotic ef- 
fects such as Majorana fermions [T31 [2] and magnetic 
monopole [U] are predicted to exist as results of the topo- 
logical SSs. 

The low-temperature scanning tunneling micro- 
scope (STM) and spectroscopy (STS) provide a direct 
way to study the SSs through probing the local density of 
states (LDOS) oscillations in the vicinity of impurities or 
step edges The quasiparticle interference (QPI) pat- 
terns induced by non-magnetic impurities on the surface 
of Bia;Sbi_a; [17] and Bi2Te3 [18], together with subse- 
quent theoretical analysis [TOti24| . demonstrated the ab- 
sence of backscattering for the topological SSs. Mean- 
while, the LDOS oscillations of SSs near step edges on 
Bi2Te3 showed a power- law decay with index —1 in a 
certain energy range |25j . compared to —1/2 for the con- 
ventional two-dimensional electron system (2DES) (16] . 
The faster decay of QPI once again indicates the sup- 
pression of backscattering in TIs. 

Despite the intensive investigation, a complete under- 
standing of QPI on the surface of TIs remain elusive par- 
tially due to the warping effect of the Dirac cone [8] [26] . 
The warping effect of the SSs results from not only the 
anisotropic surface band dispersion, but also the coupling 
between the surface and bulk bands. In this Letter, we 
present a general formalism to account for the complex 
scattering geometry. We propose the interference pat- 
terns are dominated by the extremal points on the con- 
stant energy contour (CEC) of 2D electron band. In ap- 



plying the theory to Bi2Te3 with strong warping effect, 
we show that the decay index varies from —3/2 to —1/2 
and finally to —1 as the energy increases. As for TIs with 
nearly isotropic Dirac cones, such as Bi2Se3 [27], the de- 
cay index is simply —3/2. Moreover, the theory can be 
extended to QPI induced by point defects and readily 
elucidate the missing of and the deviation of Fermi 
velocity in Ref. [18j. To confirm the predictions, we have 
performed STM study on both Bi2Te3 and Bi2Se3 films 
and found excellent agreement with theory. In particular, 
we have been able to obtain the decay index on Bi2Se3, 
whose interference pattern is usually too weak to extract 
the information. 

We start with a general 2D surface band with a single 
Fermi surface within the surface Brillouin zone (SBZ). 
Because of elastic scattering, the incoming surface wave 
with a wave vector k* must be scattered into the out- 
going one with k-'' on the same CEC. Assuming a step 
edge along the y-direction, the ky component of the wave 
vectors should be conserved in a scattering process, i.e. 
ky = = ky. The interference between the incoming 
and outgoing waves gives rise to the standing wave oscil- 
lation in the x-direction. The total LDOS is the sum of 
contributions from all these oscillations from the SSs on 
a CEC. For a given energy E, we can integrate over ky 
on the entire CEC and express the LDOS explicitly as 



5p{E,x) = 



2r/(l + |r|')elC/e'('='-'=^)"dfc, 



(1) 



where r is the reflection coefficient. Here the scattering 
wave function is of the form ^g^k^x+ikyy denotes the 
spin wave function). 

A pair of scattering states k* and k-^ lead to a stand- 
ing wave with a spatial period of 2n/{kl — fc^). Since the 
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FIG. 1: (color online) Illustration of scattering off the step 
edge along y-direction. The LDOS oscillation is dominated 
by scattering between the EP pairs on CEC (red dots), (a) 
The shape of CEC evolves from circle, hexagon to concave 
hexagon as energy increases in Bi2Se3. It has a single pair 
of EPs at (k, — k) and only one type of oscillation pattern 
appears at different bias voltages, (b) In Bi2Te3, as energy 
increases, CEC evolves from hexagon, concave hexagon to 
concave hexagram, and the pairs of EPs are first at (ki, — ki); 
then at (ki, -ki), (k2, -ka), (ka, -k2); then at (k2, -ks), 
(ks, — k2); finally at (k4, — ks) and (ks, — k4). Different types 
of oscillation pattern appear at different bias voltages. 



period is different for different value of ky, only the pairs 
whose periods are stationary with respect to small varia- 
tion in ky can make dominant contribution to the LDOS 
oscillations. We call such pair of points (with identi- 
cal ky) on CEC the extremal points (EPs) [28 . Other 
standing waves interfere with each other and cancel at 
large x. The spatial dependence of LDOS oscillations 
m Eq. ^ can be evaluated by expanding the relevant 



quantities around each EP, namely, let ky 



then kl — kl = Ak-r 



Sky 



^xo -t- J2n Afci:„(5fcy , r = J2i Vi^k'y, and 
XrnSky^. Here Ak^o is the characteristic wave 
vector depending on the geometry of CEC. To the leading 
order of ky, the LDOS varies at long distance as 



6p{E,x) ~ 



^ / 2r/(l + |r|2)ek/e^^''-'^^^dfc. 



'^IgilaXbcl cos{Akxox + (l)s)x (2) 

EPs 



and 
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is given 



where a=min(Z), 6=min(m), c=min(n), 
the initial phase of each pair of EPs. 
by /^dfc'fc'('^+''-'=+i)/'=e*'^'=-'='. The decay behavior of 
LDOS in Eq 



2^ is valid as long as x ^ ^^xQ- The 
decay index associated with a pair of EPs is given by 
{a + b + l)/c, which is solely determined by the proper- 
ties of the scattering wave function around the EPs. 

Now we turn to the topological SSs on Bi2Te3 and 
Bi2Se3 with a single Dirac cone near F point in the SBZ 
on each surface. The effective model describing such 
topological SSs reads HH ^ 



■H(k) = v{axky 



(Jy hx 



A 
2 



where h= 1, k± = ky ± ikx, v is the Dirac velocity, A is 
the warping parameter, and ai are Pauli matrices acting 
on spin space. For simplicity, here we ignore the particle- 
hole asymmetry as it affects the shape of Fermi surface 
little. The surface band dispersion is 



(^kx , ky 



± 



A2/c6 



(30), 



(4) 



where e± denotes respectively the upper and the lower 
energy bands touching at the Dirac point, and k = (fc, 9) 
with 6 as the angle between the wave vector k and k^- 
axis (F-M). The step edge is always along the close 
packed F-^ direction. Defining the characteristic energy 
e* = v^Jv/\ and length ^JX/v, in Fig. [l] we plot a set of 
CEC of the upper band in momentum space for Bi2Se3 
and Bi2Te3, respectively. In Bi2Se3 A = 128eV-A'^ and 
£* — 0.59eV, so that the CEC is nearly a circle from to 
0.42e* (0.25eV) [57]. We plot four representative CEC 
shown in Fig. [l|a). When the Fermi energy increases, 
the shape of CEC evolves from a circle {Ei = O.Sle*), 
more hexagon-like {E2 = 0.55e*), hexagon {E^ — 0.7e*) 
and to concave hexagon {E4 — 0.83e*). In a wide range 
of energy only a single pair of EPs exists at (k, — k), 
so the characteristic wave vector is always equal to 2k 
and c = 2. In Bi2Te3 the warping effect is stronger 
with A = 250eV-A3 and e* = 0.23eV [8 . As shown in 
Fig. [ijb), EPs evolve with the energy as follows: Single 
pair of EPs (ki,— ki) at Ei = 0.7e*; Multiple pairs of 
EPs (ki,-ki), (k2,-k3) and (k3,-k2j at E2 = lA6e* > 
Ec = 3i/Vll/9e* - 1.45e*; Two pairs of EPs (k2,-k3) 
and (k3,— k2) survive at E^ = 1.91e*, as the SSs along the 
F-M direction merge into the bulk conduction band; No 
EPs at all at E^ = 2.4s*, because the SSs separate from 
bulk one only in the very vicinity along T-K direction 
on the Fermi surface as observed in the ARPES experi- 
ment [5]. In this case, scattering between states around 
(k4,— ks) and (ks,— k4) will be dominant for LDOS os- 
cillations. Thus in Bi2Te3 the characteristic wave vector 
and the LDOS oscillation period critically depend on the 
bias. In most cases we have parameter c = 2 except for 
Fermi energy as high as -E4 (c = 1). 

For the incoming state with wave vector k*=(A;*,^^*) 
and energy e_|_(fc*), the inner product of two spin wave 
functions^^/ = sin6'^ -I- iA(P)3 sin(36'0 cos 6'Ve+(P). It 
is zero only when 6'' = as the spins of the time-reversal 
pairs (k, — k) are orthogonal. Thus, & = 1 in Bi2Se3; 
while in Bi2Te3, 6=1 for the pair (ki,— ki), and 6 = 
for other pairs of EPs. Assuming the step edge potential 
is V{x) = for a; < and V{x) = -V^ (Vb > 0) for 
X > 0, by matching the boundary condition at the edge 
the reflection coefficient can be obtained as 



(5) 



where ^(6**) = (e+(fc*)/fc* 
(3) A(fcO'sin(30-^)), e+(fcO 



-\{k'fsm{W'))l{e+{kf)/y~ 
= e+{y) - Vq and Qf {6') = 
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FIG. 2: (color online) LDOS oscillations due to the step edge scattering on the surface of Bi2Se3 and Bi2Te3. (a) The STM 
topograph of the Bi2Se3(lll) film, in which a step edge on the surface is observed, (b) Energy dispersion deduced from the 
standing waves at the step edge, (c)-(j) dl/dV maps and LDOS oscillations for several bias voltages obtained on the upper 
terrace in film shown in (a), (k)-(n) The LDOS oscillations on Bi2Te3 at 4 bias voltages. The LDOS is averaged over the width 
of the step and represented by red dots, and the fitting lines are in blue. 



-ef{-e^). Due to the constraint by TRS, r{6') = 
-r{-9'), and r{0' = 0) = 0. Thus a 1 for (k, -k) 
pair in Bi2Se3 and (ki,— ki) pair in Bi2Te3, and a = 
for other pairs in Bi2Te3. 

In short, the algebraical decay index is 3/2 for (k, — k) 
and (ki,— ki) pairs, 1/2 for (k2,— ka) and (ka,— k2) pairs, 
and 1 for (k4,— ks) and (ks,— kg) pairs. Therefore, the 
LDOS oscillations of the SSs in Bi2Se3 should decay as 
2^-3/2 wide range of energy (as long as i? < 0.85 eV), 
much faster than as in 2DES TO. On Bi2Te3 sur- 

faces, as the bias increases, LDOS oscillations decay first 
as x^^^^ [E < 0.33 cV), then as a combination of x^^^^ 
and a;~^/^, then as a;^^/^, and finally x~^ . 

To experimentally confirm the above predictions, we 
analyzed the interference fringes at the step edges on 
Bi2Se3 and Bi2Te3 surfaces. All experiments were car- 
ried out at 4.2 K in an ultrahigh-vacuum low temper- 
ature STM system (Unisoku) equipped with molecular 
beam epitaxy (MBE) for film growth. The stoichiomet- 
ric films of Bi2Se3 and Bi2Te3 were respectively prepared 
on graphene and Si(lll) substrates by MBE [301 E]- A 
typical STM image of Bi2Se3 film with a thickness of 50 
quintuple layers is shown in Fig. [2ja). We can clearly 
see the atomically flat morphology of the film and the 
step of the height of a quintuple layer. The steps are 
preferentially oriented along the three close packing (T- 
K) directions. The LDOS of electrons at energy eV is 
measured through the differential tunneling conductance 
dl/dV maps by STS. The Fermi velocity along f-M is 



deduced to be 4.8 x 10^ m/s by fitting the linear dis- 
persion curve [Fig. [2]^b)], in good agreement with the 
first-principles calculation and the ARPES data O [10] . 
Figs.[2];c)-[2]^j) exhibit the dl/dV maps on the upper ter- 
race by the step shown in Fig. [2|a) at various bias volt- 
ages, and the LDOS as a function of the distance x from 
the step. The Dirac point is at about 0.18 eV in STS, so 
the energy of the surface electron counted from the Dirac 
point in Fig.|2];j) is 0.38 eV (0.64e*). The best fit to the 
LDOS oscillations is given by Sp cx cos(Afc3,oX -I- (j))x~^^'^ 
as predicted. The suppression of backscattering and the 
circle-like shape of CEC lead to a much faster decay of 
LDOS in Bi2Se3 than that in 2DES 

The case of Bi2Te3 is even more interesting because of 
the strong warping effect. The LDOS oscillations at 4 
different bias voltages on Bi2Te3 film are shown in Figs. 
[2jk)-(n), and the corresponding dl/dV maps were re- 
ported in Ref. [T5] already. The Dirac point is estimated 
to be at —0.21 eV in STS fitting, thus the energies of the 
surface electron in Figs.[2]jk)-[2]jn) counted from the Dirac 
point are respectively 0.20 eV (0.87e*), 0.41 eV (1.78£*), 
0.51 eV (2.18£*), 0.61 eV (2.60e*). To extract the de- 
cay behavior from the LDOS data, for each energy we 
have compared the fittings with different decay indices. 
The best fitted LDOS oscillations clearly demonstrate 
the crossover of the decay from x~^^^ to a combination 
of and a;~^/^, then to x~^^^, and then to x~^ as the 

bias voltage increases, which agree well with our theoret- 
ical predictions. Besides, as predicted theoretically and 
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FIG. 3: (color online) The scattering geometry due to a non- 
magnetic impurity on Bi2Te3. When the CEC is a concave 
hexagram, the high DOS regions are denoted by bold lines 
along the direction of T-K. Two kinds of characteristic scat- 
tering wave vectors include the forbidden qi (dashed arrow) 
and the allowed (solid arrow) which connects a pair of EPs. 

shown in Figs.[2jl)-[2]jn), the LDOS oscillates with shorter 
period as increasing the bias. It is noticeable that due to 
the strong warping effect of the CEC, the decay rate of 
LDOS oscillations on Bi2Te3 is not always faster than in 
2DES. 

In the last part, we extend our general formalism 
for the step edge scattering to the non-magnetic im- 
purity scattering on Bi2Te3. Suppose an incident wave 
oc e''^ ''^/V^ is scattered into oc f{9)e^^ I ^ by 
the impurity potential oi U — V^biy)^ then the LDOS 
oscillates as cx f{9)e'^^^ -k') r^j, with f{0) denoting 
the scattering amplitude. The characteristic wave vector 
of the QPI pattern at large distance is obtained when 
gz(k-f-k') r jg stationary for certain direction f. With the 
concave hexagram Fermi surface of Bi2Te3 as shown in 
Fig. [3j there exist only two kinds of characteristic wave 
vectors: qi along T-K and q2 along T-M direction. Ob- 
viously, the qi connects a pair of TRS states whose scat- 
tering is forbidden. The q2 connects a pair of states at 
EP, which dominates the spatial decay. In such scatter- 
ing geometry, numerically we find q2 vary linearly with 
the energy and q2 = 1.5A; where k is the length of F-ky. 
Together with the STM data in Rec. [TE] we obtain the 
Dirac velocity along T-K as u = 4.15 x 10^ m/s , in good 
agreement with ARPES result {v = 4.05 x 10^ m/s) 8J. 
In Bi2Se3, the CEC is circle- like up to 0.25 eV and the 
characteristic scattering wave vector is always along the 
diameter. Therefore, we expect the Fourier transforma- 
tion of LDOS on Bi2Se3 surface is ring-like. 

In conclusion, our theoretical and experimental investi- 
gations indicate that the LDOS oscillation on the surface 
of TIs is generally determined by the scattering between 
surface states around the extremal points on Fermi sur- 



face, either by step edges or by non-magnetic impuri- 
ties. We have directly observed different standing wave 
patterns caused by scattering off a step on Bi2Te3 and 
Bi2Se3 with various warped surface bands, which, to- 
gether with the decay indices at different bias voltages, 
clearly demonstrate the 2D Dirac nature of topological 
surface states. 

We thank S. C. Zhang, Y. Y. Wang, and R. B. Liu for 
helpful discussion. This work is supported by the NSFC 
Grant No. 11074143, and the Program of Basic Research 
Development of China Grant No. 2011CB921901. 



* xc@mail.tsinghua.edu.cn 
t bfz@mail.tsinghua.edu.cn 
[1] X. L. Qi and S. C. Zhang, Phys. Today 63, No. 1, 33 
(2010). 

[2] J. E. Moore, Nature 464, 194 (2010). 
[3] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 
(2010). 

[4] X. L. Qi and S. C. Zhang, arXiv: 1008.2026 [Rev. Mod. 

Phys. (to be published)]. 
[5] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 

106803 (2007). 

[6] J. E. Moore and L. Balents, Phys. Rev. B 75, 121306 

(2007) . 

[7] D. Hsieh et al, Nature 452, 970 (2008). 

[8] Y. L. Chen et al, Science 325, 178 (2009). 

[9] H. Zhang et al., Nature Phys. 5, 438 (2009). 
[10] Y. Xia et al, et al., Nature Phys. 5, 398 (2009). 
[11] D. Hsieh et al, Phys. Rev. Lett. 103, 146401 (2009). 
[12] Y. Zhang et al. Nature Phys. 6, 584 (2010). 
[13] L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 

(2008) . 

[14] X. L. Qi et al, Phys. Rev. Lett. 102, 187001 (2009). 

[15] X. L. Qi et al. Science 323, 1184 (2009). 

[16] M. F. Crommie, C. P. Lutz, and D. M. Eigler, Nature 

363, 524 (1993). 
[17] P. Roushan et al. Nature 460, 1106 (2009). 
[18] T. Zhang et al, Phys. Rev. Lett. 103, 266803 (2009). 
[19] X. Zhou et al, Phys. Rev. B 80, 245317 (2009). 
[20] W. C. Lee et al, Phys. Rev. B 80, 245439 (2009). 
[21] Q. H. Wang, D. Wang, and F. C. Zhang, Phys. Rev. B 

81, 035104 (2010). 
[22] H. M. Guo and M. Franz, Phys. Rev. B 81, 041102 

(2010). 

[23] R. R. Biswas and A. V. Balatsky, arXiv: 1005.4780. 
[24] R. R. Biswas and A. V. Balatsky, Phys. Rev. B 83, 
075439 (2011). 

[25] Z. Alpichshev et al, Phys. Rev. Lett. 104, 016401 (2010). 
[26] L. Fu, Phys. Rev. Lett. 103, 266801 (2009). 
[27] K. Kuroda et al, Phys. Rev. Lett. 105, 076802 (2010). 
[28] L. M. Roth, H. J. Zeiger, and T. A. Kaplan, Phys. Rev. 

149, 519 (1966). 
[29] C. X. Liu et al, Phys. Rev. B 82, 045122 (2010). 
[30] Y. Y. Li et al. Adv. Mater. 22, 4002 (2010). 
[31] C. L. Song et al, Appl. Phys. Lett. 97, 143118 (2010). 



